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A z ~ 3 radio galaxy and its protocluster: evidence for a 
superstructure? 
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ABSTRACT 

We present spectroscopic follow-up observations of Lyman Break Galaxies (LBGs) selected 
in the field surrounding the radio galaxy MRC 0316-257 at z ~ 3.13 (0316). Robust spec- 
troscopic redshifts are determined for 20 out of 24 objects. Three of the spectroscopically 
confirmed galaxies have 3.12 < z < 3.13 indicating that these objects reside in a proto- 
cluster structure previously found around the radio galaxy. An additional 5 objects are found 
1600 km s" 1 blue-shifted with respect to the main protocluster structure. This is in addition 
to three [Om] emitters found at this redshift in a previous study. This is further evidence that 
a structure exists directly in front of the 0316 protocluster. We estimate that the foreground 
structure is responsible for half of the surface overdensity of LBGs found in the field as a 
whole. The foreground structure is associated with a strong surface density peak 1.4 Mpc to 
the North- West of the radio galaxy and a 2D Kolmogorov-Smirnov test indicates that the spa- 
tial distributions of the 0316 and foreground galaxies differ at the 3cr level. In addition, we 
compare the properties of protocluster, foreground structure and field galaxies, but we find 
no significant differences. In terms of the nature of the two structures, a merger scenario is 
a possible option. Simple merger dynamics indicates that the observed relative velocity of 
1600 km s" 1 can be reproduced if the two structures have masses of ~ 5 x 10 14 M© and have 
starting separations of around 2.5 to 3 Mpc. It is also possible that the foreground structure 
is unrelated to the 0316 protocluster in which case the two structures will not interact before 
z = 0. 

Key words: galaxies: evolution - galaxies: high-redshift - galaxies: clusters: individual - 
cosmology: observations - cosmology: early Universe 



1 INTRODUCTION 

The exact role the environment plays in the evolution of galax- 
ies has been a long standing question in astronomy. It has been 
widely observed that the properties of galaxies depend on envi- 
ronment. High density environments in the local Universe, such 
as galaxy groups and clusters, are dominated b y red early-type 
galaxies tores sled [l980l: iButcher & Oemlerll 19841) . This is in con- 
trast with the lower density environments where blue star-forming 
late-type galaxies are more frequent. Furthermore, elliptical galax- 
ies in cluste r environments are generally older than their field 
counterparts (Clemens et al. 2006; Sanchez-Blazq uez et al] 120061 : 



Ivan Dokkum & van der M arel 20071: iGobat et al.ll2008h and the red 
sequence extends down to fainter magn itudes in denser en viron- 
ments indicating it is formed earlier there (ITanaka et al.l2005U2007l 



E-mail : kuiper @ strw. leidenuniv. nl 



2008). In addition, cD galaxies, the most massive galaxies known, 
are located exclusively in galaxy clusters. 

In order to adequately explain the differences between low- 
and high-density environments it is essential to study galaxy clus- 
ters at all epochs. By doing so it may be possible to identify 
what exact physical processes constitute the more general term 
of 'environmental influence'. Unfortunately, the search for galaxy 
clusters at > 1.5 is difficult and only a handful of spectroscop- 
ically confirmed galaxy clusters at z > 1.5 with X-ray emis- 
sion are cu rrently known (Wilson et al. 2008; Papovich et al.l2010l : 
ITanaka et al.ll2010bl : lHenrv et al.ll2010l : lGobat et al.ll201ll) . 

One of the most successful methods to push the search for 
galaxy clusters bey ond z = 2 is t argeting high-z radio galaxies 
(hereafter HzRGs. lMiTev & De Breu ck 2008). With large observed 
K band luminosities, these galaxies are thoug ht to have large stel- 
lar masses of the order of 1 11 - 10 12 M© dRocca-Volmerange et al.l 
l2004l ; ISevmour et al .120071) . In the model of hierarchical galaxy for- 
mation these massive galaxies should be located in dense environ- 
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ments and are thus possible members of galaxy cluster progenitors. 
These structures are often referred to as 'protoclusters', because 
at these redshifts galaxy clusters are likely still in the proc ess of 
format ion and therefore have not yet virialised (e.g. iKuiper et al.l 
l2011ah . 

By targeting a HzRG field with a narrowband filter chosen 
such that it contains a strong emission line at the redshift of the 
radio galaxy, it is possible to select galaxies in a narrow redshift 
interval around the radio galaxy. This has resulted in evidence that 
HzRGs indeed probe overdense regions in the ear l y Universe (e.g. 
jPascarelle eTalJl 19961: Iknopp & Chambers! 1 1 9971: IPentericci etaLJ 
2000: iKurk et al.1 l2004bl| al: IVenemans et al.1 [20071 : iMatsuda et aD 
1201 ll : IKuiper et al JI201 lbh . 

Although this method is efficient in locating overdensities, it 
is limited in the number and type of galaxies that can be selected. 
The emission line most commonly used for these searches is Lyor, 
thus only galaxies with strong Lya emission are selected. As a con- 
sequence, a large number of galaxies that reside in the overdensity 
are missed altogether. This method is therefore not suitable for in- 
depth studies that attempt to obtain a complete picture of the pro- 
tocluster. 

There is a variety of methods that are aimed at selecting high-z 
galaxies. The most well-known uses the Lyman break to select UV- 
bright star-forming galaxies (Lyman Break Galaxies or LBGs) in a 
relatively broad redshift range compare d to the na r rowba nd tech- 
nique. This method was pioneered by Steidel et al.l (1 19951) and has 
been succesfully used in many studies since then. The Lyman break 
selection technique selects a much larger sample of star forming 
galaxies than methods relying on Lya narrowband data, as only 20 
per cent of LBGs at a given luminosity also qualify as Lya emitters 
taking into account current selection criteria (ISteidel et al J |2QOOL 
1201 lh . 

The Lym an break method h a s been used only sp aringly on 
HzRG fields dlntema et al] l200d : lOverzier et all l2008h . This is 
mainly because the redshift range probed by the LBG criterion is 
~ 0.3 - 0.7, significantly larger than the redshift range spanned 
by a typical protocluster (Az ~ 0.03 or Av ~ 2000 km s -1 at 
z ~ 3). One of the HzRG fields for which the LBG selection cri- 
terion has been used is MRC 0316-257 at z - 3.13. This is one 
of the best studied HzRGs at z > 2 .5 and it has been sho wn to 
host an overdensity of Lyof emitters (jVenemans et al .1120051 here- 
after V05). A study bv lMaschietto et al.1 J2008L hereafter M08) has 
found tentative evidence for a similar overdensity of [Om] emitters. 
In an attempt to obtain a complete galaxy census of the overden- 
sity, Kuip er et al.l (120 id hereafter K10) assembled photometry of 
the field in 18 bands ranging from U band to Spitzer 8 /um. K10 se- 
lected LBGs in the 40 ' 2 field around the radio galaxy and detected 
a small surface overdensity. The most massive and actively star- 
forming LBGs were found to be located near to the radio galaxy 
indicating the presence of environmental influence. However, the 
inability to distinguish protocluster galaxies from field galaxies is 
likely to diminish any real trend in the data. 

In this work we present spectroscopic follow-up of the sam- 
ple of LBGs composed in K10. By spectroscopically confirming 
the redshifts of the individual LBG candidates we can determine 
which galaxies are truly in the protocluster and which are in the 
field. This will therefore give us a better estimate of the volume 
overdensity. Also, it allows us to compare the properties of field 
and protocluster LBGs in a fully self-consistent manner. This is 
particularly important in determining whether the protocluster en- 
vironment influences the evolution of its constituent galaxies. 

The paper is structured as follows: a brief summary of the 



sample selection of K10, a description of the data and its reduc- 
tion are given in Sect.|2] Spectroscopic redshifts and the resulting 
velocity distribution are discussed in Sect.|3]and further discussion 
concerning the presence of a possible superstructure is presented 
in Sect. HI Finally, conclusions and a future outlook are presented 
in Sect. [5] Throughout this paper a standard A cold dark matter 
(ACDM) cosmology is used, with H = 71 km s _1 , Q M = 0.27 and 
Cl A = 0.73. All magnitudes are given in the AB magnitude system. 



2 SAMPLE SELECTION & DATA 

In K10 a UVR colour criterion was introduced that was designed to 
select star-forming galaxies in the redshift range 3.0 < z < 3.3. 



U-V > 1.9, 

V-R<0.51, 

£/- V> 5.07 x(V-tf) + 2.43, 

R<26. 



(1) 



A total of 52 galaxies were found to satisfy the criterion. Photo- 
metric red shifts were derived for this sample using the eazy code 
(iBrammer et al.l2008l) and broadband photometry in 18 bands rang- 
ing from U band to Spitzer 8 /um. The initial sample was then re- 
duced to 48 by applying a photometric redshift cut of 2.8 < z pho t < 
3.5. 

K10 also constructed an additional sample of 55 potential 
LBGs (pLBGs). These objects satisfy all selection criteria, except 
that they are too blue to make the U - V > 1.9 cut. However, all 
these objects are undetected in the U band used and deeper U band 
data may yield redder U-V colours. These objects are thus not 
strictly LBGs when considering the selection criterion of K10, but 
deeper data may show that they do in fact satisfy all criteria. 

Only objects with R < 25.5 were considered for follow-up 
spectroscopy, because the continuum and absorption lines of fainter 
objects are unlikely to be detected. This reduced the samples to 29 
LBGs and 27 pLBGs, respectively. The samples were subsequently 
divided in three brightness categories: objects with R < 24.5 are 
classified as 'bright', objects with 24.5 < R < 25.0 as 'intermedi- 
ate' and objects with R > 25.0 as 'faint'. To ensure the most de- 
tections, the objects in the mask were prioritised according to their 
brightness. Further restrictions were imposed by the locations of 
the individual objects as slits in the mask are not allowed to over- 
lap. The final mask contained 13 LBGs and 11 pLBGs of which 
10 are classified as 'bright', 9 as 'intermediate' and 5 as 'faint'. 
Therefore a total of 24 protocluster candidates were observed spec- 
troscopically. One of these objects is the Lya emitter #1867 from 
V05. This galaxy has been spectroscopically confirmed to be at the 
redshift of the protocluster. 

The spectroscopy was performed with the FOcal Reducer and 
low dispersion Spectrograph (FORS2) in the mask multi-object 
spectroscopy mode (MXU) at the Very Large Telescope during the 
nights of 10 and 11 December 2010. The seeing varied during the 
two nights between 0.7" and 1.2". The width of the slits in the mask 
was 1.0". The objects were observed through the "300V" grism and 
GG435 blocking filter, with a resolution of 440. The spectral range 
covered is approximately 4500 < A < 8500 A. The pixels were 
binned 2x2, which resulted in a spatial scale of 0.25" pixel -1 and 
a dispersion of 3.36 A pixel -1 . A total of 25 exposures of 1560 
seconds each were obtained. Between the individual exposures, the 
pointing of the telescope was shifted in steps of 0.25" along the slit 
to enable more accurate sky subtraction and cosmic ray removal. 
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The total integration time per object was 39,000 seconds (10.83 
hr). 

Data reduction was performed with various irai^J routines. 
The reduction included the following steps: individual frames were 
bias subtracted and flat fielded using lamp flats. Cosmic rays were 
identified and removed before the background was subtracted. The 
background subtracted two-dimensional frames were combined 
and one-dimensional spectra were extracted. Wavelength calibra- 
tion was performed using arc lamp spectra and night sky lines in 
the science frames. The uncertainty in the wavelength calibration 
is ~ 0.3 A, which corresponds to a systematic redshift uncertainty 
ofo%~ 0.0002. 



3 RESULTS 

3.1 Redshift determination 

Spectroscopic redshifts are obtained for 20 objects out of the 24 
observed. For the other 4 objects no continuum or emission lines 
are detected. Examining the data we find artefacts in two of the 
2D spectra indicating slit defects as the possible cause of these two 
non-detections. The remaining two undetected objects either have 
a low surface brightness or are very faint (R ~ 25.5) making it im- 
possible to obtain a spectroscopic redshift. Properties of the objects 
for which a spectrum was extracted are given in Table Q] 

The spectra of the objects that do allow for a redshift determi- 
nation are shown in Fig. Q] The spectra have been plotted in the 
restframe to facilitate comparison between the different objects. 
Also, the locations of the most important spectral features have 
been marked. The 2D spectra have also been included because the 
Lya break is more obviously apparent in the 2D spectra. 

The results obtained from the spectra are summarised in Ta- 
ble O Spectroscopic redshifts based on both emission and absorp- 
tion lines are listed as outflows may affect the spectroscopic redshift 
determined from Lya (e.g. Shaple v et al.l 120031) . For the spectro- 
scopic redshift we take the mean value of the redshifts obtained for 
the individual discernible absorption features. Uncertainties listed 
in Table [T] are calculated by varying the spectra according to a nor- 
mal distribution characterised by the rms noise level. The individual 
lines are subsequently fitted again. This process is repeated 1000 
times for each of the spectral features. 

Approximately half of the objects show an emission line, 
which is assumed to be Lya. This is consistent with other spec- 
troscopic studies of LBGs at z ~ 3 (e.g. Shaplev et al. 2003; 
ISteidel et al.l201ll) . The emission line is used for determining a pre- 
liminary redshift. Based on this redshift the spectrum is searched 
for consistent absorption lines. For objects that do not show an 
emission line the redshift is determined by identifying multiple in- 
terstellar absorption lines such as Siii/11260, Cii/11335 or CiV/11549 
in combination with a possible spectral break. Almost all objects 
show either a combination of an emission and an absorption line or 
multiple absorption lines, resulting in robust redshifts. 

There is one object where the redshift is possibly ambiguous. 
Object #12 has only one identifiable absorption feature and evi- 
dence for a break. Since the Lyman break is a poor redshift indi- 
cator, this makes it difficult to set an accurate redshift. The single 
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Table 1. Properties of all LBGs and pLBGs for which spectra can be ex- 
tracted. 



Object ID 


RA 


Dec. 


Type 


R 


#1 


03 


18:01.10 


-25:35:56.1 


pLBG 


25.09 


#2 


03 


18:12.30 


-25:35:42.2 


pLBG 


24.80 


#3 


03 


18:08.72 


-25:35:22.4 


LBG 


24.11 


#4 


03 


18:11.53 


-25:35:08.2 


LBG 


24.52 


#5 


03 


18:08.94 


-25:34:59.6 


LAE/LBG 


23.77 


#6 


03 


18:05.35 


-25:34:40.8 


LBG 


24.95 


#7 


03 


18:07.75 


-25:34:26.1 


pLBG 


24.37 


#8 


03 


18:18.40 


-25:34:17.1 


LBG 


24.54 


#9 


03 


18:20.28 


-25:34:02.0 


LBG 


25.06 


#10 


03 


17:58.79 


-25:33:49.9 


pLBG 


24.12 


#11 


03 


18:08.39 


-25:33:39.6 


LBG 


24.69 


#12 


03 


17:58.93 


-25:33:27.7 


pLBG 


25.19 


#13 


03 


18:04.10 


-25:33:09.1 


pLBG 


25.42 


#14 


03 


18:07.43 


-25:32:51.2 


pLBG 


24.47 


#15 


03 


18:19.41 


-25:38:16.9 


pLBG 


24.42 


#16 


03 


18:13.19 


-25:38:06.1 


pLBG 


24.83 


#17 


03 


17:58.33 


-25:37:59.0 


pLBG 


24.94 


#18 


03 


17:59.20 


-25:37:37.9 


LBG 


23.80 


#19 


03 


18:19.73 


-25:37:26.8 


LBG 


23.76 


#20 


03 


18:03.91 


-25:37:14.1 


LBG 


25.44 


#21 


03 


18:02.78 


-25:37:05.5 


pLBG 


24.99 


#22 


03 


18:13.25 


-25:36:39.7 


LBG 


23.94 



absorption feature is very strong and based on the approximate lo- 
cation of the break, the line can be identified as either the Oi/Sin 
doublet at ~ 1303 A or Cn at 1335 A. This indicates its redshift is 
either z ~ 3.11 or z ~ 3.01. To ascertain which is the more likely, 
the 2D spectrum of #12 is compared to the 2D spectra of other 
objects with strong breaks and clear Oi/Sin and Cn features. The 
former option better resembles the other 2D spectra and therefore 
we conclude that the redshift of #12 is z = 3.1 127. 



3.2 Redshift distribution 

The full distribution of spectroscopic redshifts based on the absorp- 
tion lines is shown in the left panel of Fig. [2] The redshift of the ra- 
dio galaxy is marked by an arrow. All spectroscopically confirmed 
LBGs have redshifts consistent with 2.7 < z < 3.5. The low-z inter- 
loper rate is therefore small. Assuming that the four non-detected 
LBG candidates are low-z galaxies the worst-case success rate is 
~ 83 per cent. 

The left panel of Fig. [2] also shows a clear concentration of 
galaxies near the redshift of the radio galaxy, consistent with the 
presence of a structure. To verify this we compare the observed 
Zspec distribution to the selection efficiency curve also shown in 
Fig. |2] The curve indicates that the LBG selection criterion of K10 
is most efficient in selecting objects between 2.8 < z < 3.5. A 
Kolmogorov-Smirnov (KS) test shows that there is a probability of 
0.012 that the observed spectroscopic distribution is drawn from 
the distribution defined by the selection efficiency curve. The two 
distributions therefore differ at the 2.5<x level. 

As M08 found evidence for a possible foreground structure at 
z ~ 3.1, we take a closer look at the redshift interval 3.05 < z < 
3.20. This is shown in the right panel of Fig. [2] The distribution in- 
cludes the [Om] emitters of M08 and the relevant LBGs presented 
in this work. We correct the redshifts of the M08 [Om] emitters that 
have Lya based redshifts, because the Lya line is commonly red- 
shifted with respect to the absorption lines. This redshift is due to 
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Figure 1. ID and 2D restframe spectra of individual objects for which a redshift can be determined. The ID spectra have been scaled to a common arbitrary 
flux scale. Vertical dotted lines denote the most important spectral features in this wavelength range. From left to right these features are: Lya, Nv/U240, 
SiilU260, Oi/SiiU1303, CnA1335, SirvA1392, 1402, SmvU527 and Civ^1549. The symbols next to the ID number indicate whether the object is located in the 
0316 structure (*), the foreground structure (f ) or the field ($). 
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Figure 1 - continued 



outflows and the resonant nature of the Lya line. The applied cor- 
rection is taken to be the mean difference in redshift of all galaxies 
in our sample that show both Lya emission and absorption lines. 
This correction is Az ~ 0.005 or Av ~ 350 km s" 1 , whic h is roughly 
consistent with the shift found bv lShaplev et al.l (120031) . 

The LAEs of V05 are not included in Fig. [2]because the nar- 



rowband filter used to select this sample does not allow detection of 
objects at z ~ 3.1. This is illustrated by the dotted line in Fig.|2] The 
distribution of LAEs is thus strongly concentrated near z ~ 3.13 by 
design and including it would unfairly skew the overall distribution. 
The exception to this is LAE #1518. K10 identified this object as 
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Table 2. Redshifts, equivalent width of the Lya line and discernible spectral features of all LBGs and pLBGs for which spectra can be extracted. a Only given 
for galaxies that show Lya in emission. 



Object ID 



-£spec,Lyar 



-£spec,abs 



EW , L ya W a 



Spectral features 



#1 


3 4004+ 00003 

o 19C0+O.OOO4 

3 2306+™ 4 
-Q.QQQ5 

a OS? 1 +0.0005 
-Q.QQQ5 

3 1343 + o-^ 1 

o 19^1+0:0010 

-0.0010 


3.3965+0 


0011 

M 

mt 

88H 

881? 
88?? 
8845 
88?? 
888 
8818 
88H 
881? 
881? 
8848 
8815 

0015 


17.9 + 1.1 


Break, Lya, Oi/Sin, Siiv 


#2 


3.1282+0 


14.9 ± 0.6 


Break, Lya, Oi/Sin 


#3 
#4 
#5 
#6 
#7 


3.2303+0 
3.0442+0 
3.1266+0 
3.1219+0 
3.0324+0 


43.9 ±2.8 
5.9 ±0.7 
39.3 ±1.3 
17.8 ±1.8 


Break, Lya, Sin, Oi/Sin, Siiv 

Break, Lya, Oi/Sin 

Break, Lya, Oi/Sin, C11, Civ 

Break, Lya, Sin, C11 

Break, Sin, Oi/Sin, Cn, Civ, Sin, Civ, Fen, Alii 


#8 


- 


2.9352+0 


- 


Break, Sin, Civ 


#9 


99c 7 +0.0042 
o.z,z,j/_ 00040 


3.2181+0 


12.8 ±3.3 


Break, Lya, C11 


#10 


- 


3.1121+0 


- 


Break, Oi/Sin, Cn, Civ 


#11 


- 


2.7795+0 


- 


Break, Oi/Sin, C11 


#12 


- 


3.1127+0 


- 


Break, Oi/Sin 


#13 


1AQ9+0.0017 
J,1UJ -0.0012 


3.0988+0 


9.2 ±2.0 


Break, Lya, Oi/Sin 


#14 


- 


3.1041+0 


- 


Break, Sin, Oi/Sin, Cn 


#15 


9 091 c +0.0006 
Z "* Z1J -0.0005 


2.9109+0 


7.8 ±0.8 


Break, Lya, Cn 


#16 


- 


- 


- 


- 


#17 
#18 


99C9+O.OOO6 

J, ^ Ji -0.0006 


99QC+0.0046 

2 9865+0-0 005 
y -Q.QQ07 

3 1003+0-0016 

°- luu -0.0009 


12.6 ±1.4 


Break, Lya, Sin 
Break, Sin, Oi/Sin, Cn, Civ, Fen, Alii 


#19 
#20 

#21 


0111 c +0.0007 
J,111J -0.0007 


2.5 ± 0.5 


Break, Lya, Sin, Oi/Sin, Cn, Civ, Fen, Alii 


- 


n?^^+0-0023 

9 QQQ^+0.0006 
^•^^ -0.0007 


- 


Break, Oi/Sin, Sin 


#22 


- 


- 


Break, Sin, Oi/Sin, Cn, Civ, Fen, Alii 
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Figure 2. Left panel: Distribution of all spectroscopic redshifts obtained in this work. The dashed curve shows the selection efficiency of the LBG selection 
criterion. Right panel: Redshift distribution of [Om] emitters and LBGs that are located between 3.08 < z < 3.18. Also included is LAE #1518 of V05 which 
has been identified as an LBG. The objects in the red dashed histogram are those LBGs that are associated with the 0316 radio galaxy, whereas the blue 
dashed histogram includes [Om] emitters and LBGs that are likely in a foreground structure. The solid and dotted curves indicate the transmission curves of 
the narrowband filters used to select the [Om] emitters and LAEs, respectively. The arrow marks the redshift of the radio galaxy in both of the panels. All Lya 
based redshifts are corrected for the commonly observed shift between Lya emission and absorption lines. 



an LBG and as such it is included in the right panel of Fig. [2] It is 
also included in the subsequent analysis where possible. 

The work of M08 aimed to identify [Om] emitters in the 
0316 protocluster. Spectroscopic follow-up of three of the candi- 
date [Om] emitters showed that these objects are not located at 
z ~ 3.13 as expected but at z ~ 3.1. Taking into account the correc- 
tion for the Lya redshifts, this amounts to a ~ 1700 km s _1 blueshift 
with respect to the radio galaxy. The confirmation of an additional 
5 objects at z ~ 3.1 presented in this work further strengthen the 
notion that a structure exists in front of the 0316 protocluster. 

The question is whether the structures are separate or whether 
they belong to one larger protocluster. When we look in detail at 



the LAE distribution of V05 we find that the latter is a possibility. 
The velocity dispersion of this sample is 535 km s _1 and the me- 
dian redshift is similar to the median redshift of the red distribution 
shown in Fig. [2] If we use a KS test to compare this to the expected 
distribution based on the transmission of the Lya narrowband filter 
(dotted line), then we find a probability of 1.2 x 10" 6 that the LAE 
distribution is drawn from the expected distribution. This indicates 
that the blueshift of the LAE distribution with respect to the central 
redshift targeted by the narrowband filter is real. This in turn im- 
plies that the distribution of LAEs may extend to lower redshifts, 
but that these objects have been missed due to the location of the 
narrowband filter. 
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On the other hand, a normalized tail index (IBird & Beer si 

1 19931) of 0.66 implies that the composite distribution more closely 
resembles a uniform distribution rather than a single Gaussian. For 
this work we will follow the approach of M08 and assume that the 
foreground structure is a separate structure. More evidence for this 
will be provided in Sect. 14.21 

Assuming that there are two subgroups we find mean redshifts 
of z = 3.1039 and z = 3.1262 which implies a velocity difference of 
~ 1620 km s _1 . Using a Gapper scale estimator (B eers et al.lll990h 
we find a velocity dispersion of 965 ±112 km s _1 for the composite 
distribution and individual velocity dispersions of 492 ±120 km s _1 
and 364 ± 120 km s _1 for the blue- and redshifted subgroups, re- 
spectively. 

Although we cannot determine, based on the redshift distribu- 
tion, whether the two groups are separate structures or one larger 
protocluster, for the remainder of this work we will refer to the 
objects associated with the possible foreground structure as fore- 
ground objects, whereas those objects at z ~ 3.13 will be clas- 
sified as 0316 galaxies. Objects not associated with either of the 
z~ 3.1 -3.13 structures will be referred to as field galaxies. 




Aa (arcmin) 



4 DISCUSSION 

4.1 A possible superstructure and implications for the 
overdensity 

In K10, the LBG surface density of the 0316 field was found to be 
a factor 1.6 ± 0.3 larger than the control field used. Using a volume 
argument and assuming that the surface overdensity is caused by a 
single structure connected to the z = 3.13 radio galaxy, K10 trans- 
lated this surface density to a volume density which is 8 ± 4 larger 
than the field density. 

Based on the surface overdensity, one expects spectroscopic 
follow-up to reveal approximately 1 out of 3 objects to be in the 
0316 protocluster. Instead, 3 out of 20 objects are found to be asso- 
ciated with the radio galaxy and an additional 5 objects are found to 
be located at ~ 3.1. We therefore confirm the surface overdensity, 
but it is not solely caused by the 0316 protocluster. The volume 
overdensity around the 0316 radio galaxy found by K10 is thus not 
as large as previously estimated and must be adjusted. 

To correct the volume overdensity we must assess what frac- 
tion of the overdensity is due to the foreground structure. Due to the 
small samples that are considered, this cannot be more than a rough 
estimate. In the sample of 13 [Om] emitters found by M08 a total of 
eight have spectroscopic redshifts; five due to overlap with the Lya 
emitter sample of V05 and three through spectroscopic confirma- 
tion of the [Om] line. The latter three were all found to be at z ~ 3. 1 . 
Taking into account that the five unconfirmed [Om] emitters may be 
in either of the two structures, the fraction of foreground objects is 
thus 25-60 per cent. 

In this work, five out of eight objects (or ~ 60 per cent) are 
found to be in the foreground structure. Based on these numbers we 
assume that half of the surface overdensity can be attributed to the 
foreground structure. We therefore estimate the volume overdensity 
of the 0316 protocluster to be ~ 4, rather than 8. This is very similar 
to the overdensity of Lya emitters of 3.3^ found in V05. Based 
on these numbers we also expect that the foreground structure is 
similar in richness and mass as the 0316 protocluster. 



Figure 3. Spatial distribution of all spectroscopically confirmed objects that 
reside in the 0316 protocluster (red symbols) or the foreground structure 
(blue symbols). Lya emitters are denoted by plus signs, [Om] emitters by 
triangles and spectroscopically confirmed LBGs/pLBGs by diamonds. Also 
shown are the locations of the unconfirmed LBGs and pLBGs as small open 
circles. The location of the radio galaxy is marked by the star. The dotted 
lines denote the size and location of the narrowband image used for the 
detection of the [Om] emitters. The contours indicate the surface density of 
unconfirmed LBGs and the LBGs that are shown to be either at z ~ 3.10 or 
z~3.13. 



4.2 Spatial distribution 

Figure [3] shows the spatial distribution of all objects that are spec- 
troscopically confirmed to be either in the 0316 protocluster or in 
the z = 3.1 foreground structure. This includes the 32 Lya emitters 
of V05 and the three [Om] emitters of M08. The blue objects are 
those identified to be in the foreground structure. 

Interestingly, four out of five of the foreground LBGs are lo- 
cated in a small area in the North- West region of the field. This 
specific region is also mostly devoid of z = 3.13 Lya emitters. The 
North- West region is thus dominated by foreground objects. Eight 
of the non-confirmed LBGs are also located in that general region, 
four of which are strongly clustered. 

The distribution of the foreground [Om] emitters does not re- 
flect this apparent concentration of objects. This is, however, due 
to the small size of the narrowband image used for the detection 
of these objects, as illustrated by the outline shown in Fig. [3] This 
limits the location of the foreground [Om] emitters to the central 
region of the field. 

To further illustrate the subclustering, Fig.|3]also shows LBG 
surface density contours. The LBGs that have been shown to be 
field galaxies are not considered, but LBGs that have not been spec- 
troscopically confirmed or have been confirmed to be in either of 
the two z ~ 3 structures are included. The contours have been ob- 
tained using a grid with a gridsize of 3". For each cell of the grid 
the surface density was calculated by determining the distance to 
the Mh nearest neighbour and subsequently using cr - N/nrjj. The 
contours shown have been obtained with N = 8. The resulting sur- 
face density map is smoothed with a smoothing length of 0.5'. A 
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clear peak in the surface density map is located in the North- West 
region, near the concentration of foreground objects and ~ 3' or 
~ 1.4 Mpc from the radio galaxy. A second, less pronounced peak 
in the surface density is located in the centre of the field, near the 
radio galaxy. This provides further evidence that the foreground 
structure is offset from the 0316 protocluster. 

The significance of the spatial subclustering can be quantified 
using a 2D KS test. First we determine whether the distribution of 
LBGs in the foreground structure is consistent with being drawn 
from a random distribution. We find a probability of 0.033, which 
implies that the distribution is different from random at the ~ 2cr 
level. 

When the spatial distribution of foreground objects is com- 
pared to that of the 0316 objects, we find a probability of 0.0034 
that both originate from the same parent distribution. The distribu- 
tions of the two structures therefore differ at the ~ 3<x significance 
level. The foreground [Om] emitters have not been taken into ac- 
count in this comparison. The foreground structure thus seems to 
be offset from the 0316 structure. We consider this evidence that 
the two structures are two separate groups of galaxies and not one 
single protocluster. 



4.3 Influence of environment on galaxy properties 

The influence of environment on galaxy evolution is an important 
topic in present day astronomy. Local galaxies in dense environ- 
ments are generally older, redder and have lower star formation 
rates (SFRs) than those in less dense environments. There has, 
however, been mounting evidence that the decrease in star forma- 
tion obser ved locally in dense regions turns around at earlier cos- 
mic times Jilbaz et al.ll200l ICooper et al.ll2008l : iTran et alJEoiol : 
Iffilton et al.ll2010l : IPopesso et al.ll201ll) . 

Protocluster fields make excellent targets for studying envi- 
ronmental effects at z > 2 and several studies have presented ample 
evidence tha t enviro nment influences galaxy properties at z ~ 2. 
iTanaka etal] (l2010al) showed that galaxies in the z ~ 2.15 pro- 
tocluster around PKS 1138-262 assembled their mass earlier than 
field galaxies. lHatch et al. (201 1) found that Ha emitters in the pro- 
tocluster around the radio galaxy 4C+ 10.48 at z = 2.35 are twice as 
massive as their field counterparts. Similarly, Stei del et al] d2005h 
showed that galaxies in a serendipitously discovered protocluster 
at z = 2.3 are approximately twice as old and twice as massive as 
their field counterparts. 

For the 0316 protocluster, no significant differences have been 
found between the field and the protocluster galaxies in terms of 
mass or SFR (K10). There are, however, trends of decreasing mass 
and SFR with increasing distance from the radio galaxy. One of the 
main problems is that the field interlopers in the LBG sample of 
K10 possibly dilute any differences that may be apparent in a pure 
protocluster sample. With the spectroscopy presented in this work 
a first division between protocluster LBGs and field LBGs can be 
made. 




4.3.1 SED fitting 

We first examine the properties of the LBGs in different envi- 
ronments by fitting their spectral energy distributions using pho- 
tomet ry presented in K 10. For this we use th e fast SED fitting 
code ( Kriek et al.l2009h in combination with the lBruzual & Charlod 
(2003) evolutionary population synthesis models. Originally, all 
objects were assumed to be located at z = 3.13. The addition of 
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Figure 4. Stellar mass distributions of galaxies residing in the field (black), 
the 0316 protocluster (red dashed) and the foreground structure (blue 
dashed). 



spectroscopic redshifts allows for a fully self-consistent compari- 
son between field and protocluster galaxies. The free parameters in 
the fitting routine are the age, mass, SFH and the extinction by dust, 
but as in K10 we only focus on the stellar mass, because this is the 
only property that can be determined with reasonable accuracy. 

Briefly summarising the details of the SED fitting process em- 
ployed in K10: we consider exponentially declining SFHs with de- 
cay times, r, ranging from 10 Myr to 10 Gyr with steps of 0.1 dex. 
The ages we consider range from log(age/yr)=7 to the age of the 
Universe at z ~ 3.13 which is log(age/yr)=9.3. The ICalzetti et al J 
(2000) extinction law is used for the internal dust extinction, with 
Ay ranging from to 3 with steps of 0.1. For all cases a Salpeter 
mass function and solar metallicity are assumed. 

The SED fitting results show a marginally larger mean stel- 
lar mass of 4.8 x 10 10 M© for the 0316 galaxies compared to 
2.9 x 10 10 M Q for the foreground structure and 1.7 x 10 10 M for 
the field galaxies. When we combine the 0316 and foreground sam- 
ples we obtain a mean mass of 3.7 x 10 10 M© for LBGs in dense 
environments. However, the small samples considered in this work 
imply that this difference is not significant. This is also apparent 
from the stellar mass distributions shown in Fig. [4] The small sam- 
ples make it impossible to distinguish the distributions. KS tests 
also reflect this, yielding probablities of 0.75 - 0.9 that the various 
distributions are drawn from the same parent distribution. There is 
therefore no discernible difference in stellar mass between the field 
and protocluster populations. 



4.3.2 Stacked spectra 

The presence or absence of environmental dependence is studied 
further using the stacked spectra. In Fig. [5] we show a series of 
stacked spectra of the galaxies in each of the environmental cat- 
egories. These stacked spectra have been obtained by shifting the 
individual spectra to a common restframe wavelength scale using 
the absorption line redshifts. The individual spectra are then scaled 
to the same mean flux level in the restframe wavelength range 
1300 < /l rest < 1500 A and subsequently added together. Here the 
observed wavelength range 5565 < A < 5590 A is excluded due to 
the presence of strong night- skyline residuals. Since we are dealing 
with small samples no other outliers are excluded in the stacking 
process. 
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Figure 5. Stacked LBG spectra including, from top to bottom: 0316 objects 
with 3.12 < z < 3.13, foreground structure objects with 3.1 < z < 3.12, 
combined sample of 0316 and foreground objects, field objects and all spec- 
troscopically confirmed objects. Dotted vertical lines indicate the location 
of the most important spectral features as in Fig.[T] Spectral features above 
1600 A are Fen/11608, Hen/11640 and AlnA1670. 



The properties of all detectable spectral lines in the stacked 
spectra are listed in Table |3] Uncertainties on the derived proper- 
ties were obtained by repeating the stacking process, but with a 
number of the spectra replaced by randomly drawn spectra from 
the same sample. This is to obtain a measure of the intrinsic scat- 
ter between the different spectra. For the field galaxies we replace 



three of the spectra, for the combined 0316+foreground sample we 
replace two spectra, whereas for the 0316 and foreground struc- 
tures only one spectrum is replaced. In total 10, 15, 20 and 50 fake 
spectra are constructed for the 0316, foreground, 0316+foreground 
and field galaxies, respectively. The fake stacked spectra are subse- 
quently varied according to their rms noise after which all proper- 
ties are recalculated. For each fake spectrum this process is repeated 
100 times. The standard deviations of the subsequent distributions 
are taken as \cr uncertainties. The Lya FWHM are measured from 
stacked spectra using the Lya redshifts rather than the absorption 
line redshifts. Using the latter increases the FWHM by a factor of 
1.5-2. 

The top panel of Fig. \5\ shows the stacked spectrum of the 
three objects identified to be in the 0316 protocluster. The main 
feature of the spectrum is the strong Lya emission with EW = 
26 .4 + 3.8 A, but this is mostly driven by the LAE included in the 
sample. Removing this galaxy lowers the E W to 17.9 A. The field 
galaxies, in the fourth panel, also show Lya emission, but with an 
EWq - 7.4 + 2.1 A it is on average not as strong as in the 0316 
protocluster. 

The average spectrum of the foreground galaxies shows no 
Lyof emission, but considering the small sample size we are un- 
able to determine w hether this is a significant difference. In fact, 
IShaplev et alj (2003) show that approximately half of all LBGs 
show Lya in emission and the other half shows Lya in absorption. 
With sample sizes of 3 for the 0316 protocluster and 5 for the fore- 
ground structure, the chance that the observed difference between 
the 0316 and foreground structure is a mere statistical fluctuation 
should be considered significant. We can therefore draw no con- 
clusions based on this sample, but It would be interesting to see 
whether this difference persists in a larger sample. 

We can slightly alleviate the problem with the sample size by 
combining the 0316 and foreground structure samples and compar- 
ing this to the field population. Any difference would be a strong 
indication of environmental influence at z ~ 3. The stacked spec- 
trum is shown in the third panel of Fig. [5] and the relevant proper- 
ties are listed in Table [3] The spectrum shows Lya emission with 
EWq - 11.3 + 3.7 A. This is consistent with the field population 
within the l<x uncertainties. The strength of the absorption lines 
of the combined sample are also consistent with that found for the 
field population. The composite spectrum of galaxies in the z ~ 3 
overdense structures therefore does not differ from field galaxies, 
so we find no evidence for environmental influence at z ~ 3. 

The properties listed in Table |3] can be compared to the re- 
sults of IShaplev etaD J2003L hereafter S03). In S03 the spectra of 
~ 1000 LBGs were stacked to perform a detailed study of the av- 
erage properties of these galaxies. In ge neral, the LBG properties 
in Table [3] are similar to those in Shaple v et al.1 (12003b . There is a 
velocity difference of ~ 400 - 900 km s _1 between the Lya line and 
the absorption lines, where the Lyor line is redshifted with respect 
to the absorption. This is also seen in the individual LBG spectra 
of S03 and is indicative of outflows. The LBGs in the 0316 field 
therefore also show evide nce of outflows. 

IShaplev et al.l (2003) divided their sample of LBGs into four 
bins based on Lya equivalent width, ranging from Lya in ab- 
sorption to strong Lya emission. LBGs with strong Lya emission 
were found to have weaker low-ionisation lines, bluer UV slopes 
and smaller kinematic offsets between Lya and interstellar absorp- 
tion lines. Based on the stacked spectra presented here we can 
make a similar division between the various samples. For this pur- 
pose we will only consider the field population and the combined 
0316+foreground population as the individual 0316 and foreground 
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Table 3. Properties of the emission and absorption lines found in the stacked spectra of the 0316 protocluster, the foreground structure and the field. Restframe 
equivalent widths are taken to be positive for emission lines and negative for absorption lines. Velocity offsets are given with respect to the Oi/Sin doublet. 
FWHM values are corrected for the instrumental resolution. a Based on a stack of the three objects that do show Lya emission. b These values are obtained 
from stacked spectra created using the Lya redshift where available. c This value cannot be constrained and is therefore not listed. d The UV slope ft is 
calculated using the R and / band data of K10. 





0316 


Foreground 


0316 +Foreground 


Field 


Av Lya (kills X ) 

EW , L ya (A) 

FWHM Lyff (kms~y 


+442 ± 142 
26.4 ± 3.8 
561 ±118 


+734 ± 225 a 
-13.3 ±5.3 

_c 


+451 ± 156 
11.3 ±3.7 
493 ± 122 


+396 ± 162 

7.4 ±2.1 
803 ± 241 


AvsmCkms" 1 ) 
EW om (A) 
FWHMsmCkms- 1 ) 


- 


+308 ± 178 
-1.9 ±0.8 
400 ± 352 


+340 ± 167 
-1.5 ±0.8 

c 


+281 ± 152 
-2.2 ± 0.8 
558 ±212 


Avoi/sniCkms" 1 ) 

£Wb,oi/sin (A) 
FWHMoi/smCkms- 1 ) 



-2.1 ±0.7 
365 ± 233 



-4.5 ±1.1 

799 ± 233 



-3.5 ±0.9 

708 ± 223 



-2.4 ± 0.5 
623 ± 201 


A vcii (kms -1 ) 
EW ,cii (A) 
FWHMcuCkms- 1 ) 


-263 ± 169 
-3.0 ±1.1 
726 ± 309 


+ 103 ±360 
-3.1 ±1.1 
1802 ±545 


-117 ±202 
-2.7 ± 0.8 
1127 ±559 


-61 ± 188 
-1.8 ±0.6 
509 ± 205 


Av S iiv(kms _1 ) 
EW ,snv (A) 
FWHMsnv (km s" 1 ) 


- 


- 


+ 151 ±196 
-1.4 ±0.8 

_c 


-61 ± 204 
-0.9 ± 0.5 
102 ± 226 


AvsniCkms" 1 ) 
EW ,siii (A) 
FWHMsmCkms- 1 ) 


- 


-75 ± 263 
-2.1 ±0.9 
711 ±266 


- 


+88 + 216 
-1.7 ±0.5 
315 ±231 


Avciv (km s _1 ) 
EWo,civ (A) 
FWHM C iv(kms- 1 ) 


-369 ± 267 

-5.3 ±1.5 
1881 ±414 


-194 ±220 
-3.1 ±0.9 
800 ± 334 


-213 ± 226 
-3.4 ±1.1 
1408 ± 607 


-120 ±324 

-2.7 ± 0.8 
887 ± 286 


Av F eii(kms _1 ) 
EW , ¥ell (A) 
FWHMpeiHkms- 1 ) 


- 


- 


- 


+24 ± 252 
-1.2 ±0.5 
414 ± 278 


pd 


-1.7 ±0.2 


-0.8 ± 0.4 


-1.2 ±0.3 


-1.3 ±0.1 



samples are too small to make a meaningful comparison. Both pop- 
ulations fall into the moderate emission category or group 3 of S03. 

Little difference is seen between the composite field spectrum 
and the results of S03. All EW values are consistent within \cr 
with the properties of group 3 in S03. Since we make no distinction 
in EW Lya when stacking the spectra we also expect this spectrum to 
match closely to the full LBG stack of S03. Indeed, all absorption 
line equivalent widths are fully consistent with the average LBG of 
S03. 

When comparing the 0316+foreground objects with group 3 
of S03 we see that the Oi/Sin doublet and the Cn line of the com- 
bined sample is stronger. This is partially due to the inclusion of 
#12 in the stack, which has an exceptionally strong Oi/Sin dou- 
blet. Removing this object from the summed spectrum reduces the 
equivalent width to -2.7 + 0.6 A which is formally consistent with 
the results of S03. This, however, does not explain the strong Cn 
feature. The other spectral features are consistent with SO 3. 

Following the physical picture presented by S03, stronger ab- 
sorption lines may be explained by a larger covering fraction of 
the outflowing gas. This would, however, also diminish the Lya 
flux. Since the 0316+foreground sample does show significant Lya 
emission there must be something compensating for the larger cov- 
ering fraction. This could be related to a lower than expected dust 
content. The UV slopes of galaxies are sensitive to dust content, 



but the values listed in [3] show no significant difference between 
this work and S03. Furthermore, we must consider that the samples 
used in this study are much smaller than the samples presented in 
S03. This could indicate that the strong Cn absorption line is due to 
statistical fluctuation. It is therefore necessary to increase the num- 
ber of spectroscopically confirmed galaxies in order to put proper 
constraints on any possible differences between the various popu- 
lations in this work and the results of S03. 

The stacked spectra can also be used to determine whether 
IGM absorption blueward of the Lya line is more prevalent in either 
of the structures or in the field. To do this we assess the mean flux 
level for /l rest < 1185 A and compare it to the mean flux level for 
/l rest > 1280 A. The ratio of these flux levels is highest for 0316 at 
0.67 + 0.02, whereas the foreground and field galaxies show ratios 
of 0.53 + 0.02 and 0.50 + 0.01, respectively. The combined sample 
of the two structures yields 0.59 + 0.01. The Lya break is thus less 
pronounced in the protocluster galaxies indicating that there is less 
IGM absorption in the overdense structures. 



4.4 Interacting or unrelated structures? 

The presence of a foreground structure is not the first indica- 
tion that HzRG- sele cted protoclusters are part of superstructures. 
iKuiper et al] (1201 lah has shown that the well- studied protocluster 
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around PKS 1138-262 (1138) atz = 2.15 exhibits a broad bimodal 
velocity structure which has been independently found for both the 
megaparsec scale structure and the central kiloparsec scale struc- 
ture. This is best explained by a line-of- sight merger scenario of 
two massive halos. Could this also be the case for the 0316 proto- 
cluster and its foreground companion? 

In Sect. 14.21 we show that the spatial distributions of the two 
structures in the 0316 field are not drawn from the same distribution 
and Fig.[3]indicates a projected separation of ~ 1.4 Mpc. If there is 
a merger it is not along the line of sight. 

The velocity difference between the 0316 protocluster and the 
foreground structure is ~ 1600 km s _1 , but this is only the line- 
of- sight velocity component. If this is indeed a merging or in- 
teracting system, then additional transverse velocity components 
may be present. This implies that the true relative velocity may 
be larger than 1600 km s _1 . The relative velocity is thus simi- 
lar to t he 1600 km s _1 found in the 1138 system. Kuiper et al. 
(2011 a|) showed that the Mille nnium simulation ( Springel et al. 
120051 ; iDe Lucia & Blaizotl 120071) could reproduce such a velocity 
difference, but only f or the largest halo m asses. Doing the same 
analysis at z ~ 3 as in Kuiper et al.l (l2011al) reveals no such merg- 
ers in the Millennium simulation. 

In order to determine whether such a merger is possible at z ~ 
3, we calculate how the relative velocity evolves with decreasing 
distance d in the case of tw o merging mas sive halos. For this we use 
the equations described in|Sarazin|(2002). Conservation of energy 
dictates 
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with m = M\M 2 I{M\ + M 2 ) and d the separation between the 
structures when they drop out of the Hubble flow. For simplicity it 
is assumed that the transverse velocity is zero. This yields 



v= A /2G(M 1 + M 2 )(i--j- 



(3) 



In Fig. [6] we show how the velocity increases with decreasing 
distance in a merger scenario. The curves shown are for a variety of 
values for d . The largest value for d was chosen such that the time 
it takes to reach a relative velocity of 1600 km s _1 is equal to the 
age of the Universe at z = 3.13. We also consider two specific halo 
masses, but in all cases it is assumed that the 0316 and foreground 
structure are of equal mass. We see that for masses of 10 14 M© a 
velocity of ~ 1600 km s _1 is only reached at small separations of 
the order of < 0.6 Mpc. For masses closer to the estimated mass of 
the 0316 protocluster, the distance at which 1600 km s _1 is reached 
ranges between 1.0 < d < 1.6 Mpc. 

As we have shown in Sect. 14.21 the projected distance between 
the structures is of the order of ~ 1.4 Mpc. If the two structures 
are interacting then the true distance is likely larger. The relative 
velocity of 1600 km s _1 is also likely a lower limit due to projection 
effects. If the merger scenario is possible we therefore expect the 
curves to cross through the upper right quadrant of Fig. [6] There 
are two curves that meet this requirement, with starting separations 
of 2.5 and 2.95 Mpc and masses of 5 x 10 14 M Q . It is therefore 
possible that the system is undergoing a merger, but only if the two 
structures are massive and the starting separation is roughly 2.5 to 
3 Mpc. 

The other option is that the 0316 protocluster and the fore- 
ground structure are 'unrelated' structures in the Hubble flow 
at the time of observing. Assuming that the Hubble flow domi- 
nates the relative motions of the two structures we use the Hub- 
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Figure 6. Evolution of the relative velocity of two massive structures with 
distance for a variety of starting distances do and halo masses. The starting 
separations are 1.5, 2.0, 2.5 and 2.95 Mpc. These are denoted by the solid, 
dotted, dashed and dash-dotted curves respectively. For each starting dis- 
tance two halo masses are considered. The halo mass ratio is fixed at 1:1 
and the massed considered are 1 and 5 x 10 14 M© with the steeper curves 
corresponding to the larger halo masses. The horizontal and vertical lines 
indicate the relative velocity of 1600 km s _1 and the projected separation 
between the two structures as observed in the 0316 field. 



ble law to find that a velocity difference along the line of sight of 
1620 k m s -1 implies a dista nce of at least 23 comoving Mpc. Fol- 
lowing iBahcaHetal] J2004I) we assume that the mean distance be- 
tween two clusters in the Hubble flow at z ~ 3 is ~ 50 Mpc. What 
is the chance of encountering such a line-of-sight alignment as pos- 
sibly witnessed here? In a field of 100 x 100 Mpc 2 we expect a total 
of ~ 15 structures. The 0316 field covers approximately 9 Mpc 2 . So 
the probability of having one additional structure directly in front 
of the main structure is ~ 1 per cent. If we lower the mean distance 
to ~ 30 Mpc, then the chance increases to ~ 4 per cent. These 
probabilities are small, but they are not negligible. It is thus not 
unreasonable to find an unrelated foreground structure and the pos- 
sibility that this is a change alignment of two unrelated structures 
in the Hubble flow cannot be rejected. 

If the two structures are a chance alignment then they are sepa- 
rated by ~ 23 Mpc. This does not preclude that at some later epoch 
the two structures will interact. We try to determine whether two 
massive haloes separated by ~ 23 Mpc at z ~ 3 will merge before 
z = by using the Millennium simulation. Only the most mas- 
sive haloes at z ~ 3 are considered, because the estimated mass of 
the 0316 protocluster is large (~ 5 x 10 14 M Q , V05) and the fore- 
ground structure seems to have a similar mass. Using a lower limit 
of 2x 10 13 M© we find a total of 20 halo pairs with relative distances 
between 20 and 25 Mpc. For the majority of the pairs the separation 
decreases, but none reach separations smaller than 8 Mpc. None of 
the pairs therefore interact or merge before z = 0. Thus, if the two 
structures in the 0316 field are not interacting at z ~ 3, then it is 
likely that there will be no interaction between the structures on the 
timescale of ~ 10 Gyr. 



5 CONCLUSIONS 

We have presented spectroscopic follow-up of LBGs in the field 
surrounding the 0136 protocluster at z ~ 3.13. We observed a to- 
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tal of 24 LBG candidates. Using spectroscopic redshifts we distin- 
guish between field and protocluster galaxies. This in turn allows 
us to make a self-consistent comparison between the field and pro- 
tocluster galaxy samples. 

(i) We determine redshifts for 20 out of 24 objects, finding that 
all objects are located between 2.7 < z < 3.5. This implies an in- 
terloper fraction of at most ~ 17 per cent. Out of the 20 confirmed 
objects, 5 are located at z ~ 3.10 and 3 at z ~ 3.13. The number 
of 0316 protocluster objects is too small to account for the surface 
overdensity presented in K10, but is consistent with the presence of 
two structures: the 0316 protocluster at z ~ 3.13 and a foreground 
structure at z ~ 3.1. The pres ence of such a foreground structure 
was already hypothesised by iMaschietto et al] ((2008). Based on 
the redshift distribution, however, it is not clear whether the two 
structures should be considered as separate or as part of one larger 
structure. 

(ii) The spatial distribution of the foreground and 0316 LBGs 
shows two distinct density peaks: one centred on the radio galaxy 
and a stronger peak located in the North- West corner of the field. 
This latter stronger peak coincides with a concentration of fore- 
ground objects indicating that the foreground structure is not di- 
rectly in front of the 0316 protocluster. A 2D Kolmogorov-Smirnov 
test confirms this, indicating that the spatial distributions of the 
0316 and foreground LBGs differ at the 3<x level. This implies 
that the two structures are likely indeed separated and not part of a 
larger protocluster. 

(iii) The presence of the foreground structure implies that the 
volume overdensity of LBGs of K 10 is overestimated. Instead of 
the previously determined overdensity of 8, the volume density is 
only a factor ~ 4 larger than the field. We also estimate that the 
foreground structure is of similar mass and richness to the 0316 
protocluster. 

(iv) There is no systematic difference in mass between the proto- 
cluster galaxies and field galaxies. Stacking the spectra shows that 
the galaxies associated with the 0316 protocluster have stronger 
hya emission than the field galaxies, whereas the galaxies in the 
foreground structure show very little hya emission. However, con- 
sidering the limited sample size this may be due to statistical fluc- 
tuation. Combining the galaxies in the two structures in one com- 
posite sample shows that galaxies in dense environments do not 
differ from field galaxies. This implies that there is no discernible 
evidence for environmental effects on galaxy evolution at z ~ 3. 
However, the Lyman break is less pronounced in the combined 
0316+foreground sample indicating that there is less absorption by 
the IGM in these structures with respect to the field galaxies. 

(v) The relative velocity of the 0316 and foreground structures 
can be reproduced if the structures are merging and have large 
masses of 5 x 10 14 M . Alternatively, the structures may not be 
merging in which case the relative velocity translates into a radial 
separation 23 Mpc. Such a large separation at z ~ 3 means it is 
unlikely that the two structures will interact by the present day. 

The results presented here for the different samples of galax- 
ies should be considered preliminary. Further spectroscopic obser- 
vations are necessary to get a better census of which galaxies are 
in which of the structures. This will also give better constraints on 
the potential differences between the protocluster and field galaxies 
at z ~ 3. Furthermore, extra data will provide stronger constraints 
on the merger scenario considered in this work. If additional spec- 
troscopic redshifts result in smaller separations on the sky and in 
redshift space, then this may increase the likelihood of the merger 
scenario. 
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